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Abstract 

There exist quantum algorithms that are more efficient than their classical counterparts; 
such algorithms were invented by Shor in 1994 and then Grover in 1996. A lack of inven- 
tion since Grover's algorithm has been commonly attributed to the non-intuitive nature of 
quantum algorithms to the classically trained person. Thus, the idea of using computers 
to automatically generate quantum algorithms based on an evolutionary model emerged. 
A limitation of this approach is that quantum computers do not yet exist and quantum 
simulation on a classical machine has an exponential order overhead. Nevertheless, early 
research into evolving quantum algorithms has shown promise. 

This paper provides an introduction into quantum and evolutionary algorithms for the 
computer scientist not familiar with these fields. The exciting field of using evolutionary 
algorithms to evolve quantum algorithms is then reviewed. 

1 Introduction and Overview 

Quantum algorithms and evolutionary algorithms are two increasingly popular research 
fields. Although still relatively new areas of research, there is a large number of publications 
in both areas. The idea of using evolutionary algorithms to produce quantum algorithms, 
known as evolving quantum algorithms, has not been pursued to the same extent as research 
in the two fields independently. This new field of evolving quantum algorithms is in an in- 
teresting development stage with the potential to dramatically change the area of quantum 
computing. 

This is a review paper in the field of evolving quantum algorithms using genetic program- 
ming. The objectives of this paper are to 

• provide a detailed introduction to the field of quantum computing, 

• introduce the concept of evolutionary algorithms, specifically genetic programming, 
and 

• discuss the application of genetic programming to quantum algorithms, including a 
comprehensive summary of work in this new field. 

These objectives will be addressed in order; this paper is directed at readers with general 
knowledge of computer science, but not necessarily of quantum computing or evolutionary 
algorithms. 



2 Introduction to Quantum Computing 

Since Feynman's famous speech delivered in December 1959 in which he mentioned the pos- 
sibility that "sub-microscopic" computers could be built [20], the world has seen great ad- 
vances in computing power and computer miniaturization, especially in the last two decade^]. 

1 A more detailed and structured account of Feynman's idea was given in his book 1211 . 



1 



Nonetheless, there are limits to computer miniaturisation with current microchip technol- 
ogy [50] and we still, as yet, have not realised Feynman's sub-microscopic prediction. The 
desire for "sub-microscopic" computers is fueled by the appeal of the ability to significantly 
increases the efficiency of storing, copying, transmitting and processing information using a 
computer that will not occupy significant physical space. 

The two major approaches that have been proposed for sub-microscopic computing are 
quantum computing and DNA computing, which have been widely discussed theoretically 
and are the subject of a large number of theoretical and empirical studies. DNA or molec- 
ular computing, not to be confused with evolutionary or genetic algorithms, essentially 
involves using DNA molecules instead of microchips. This allows the information-processing 
capabilities of organic molecules to replace digital switching primitives and achieve this sub- 
microscopic size [50]. As this section of this paper only introduces quantum computing, more 
information on DNA computing is contained in the excellent book by Paun, Rozenberg and 
Salomaa [51] or the shorter article by Grofi [28] , 

2.1 Brief History of Quantum Computing 

Quantum computing is derived from a link between quantum mechanics, computer science 
and classical information theory 39 . Essentially quantum quantum computing can be 
viewed as developing algorithms that will run on quantum computers, which are based on 
quantum mechanics. In the 1920s, the theory of quantum mechanics was proposed: the well- 
known major contributors were Born, Dirac, Heisenberg, and Schrodinger [39]. However, 
it wasn't until 1980 that the area of quantum computing truly began, when Benioff [H] 
presented his quantum Turing machine. It proved that quantum systems could coherently 
perform computation to solve problems. 

Two years later in 1982, Feynman [T5] observed that quantum computing could not 
always be simulated efficiently on a classical computer, which lead to speculation that quan- 
tum computing could be more efficient than classical computing in some cases. In 1985, 
Deutsch [JD further developed Benioff 's quantum Turing machine and suggested that quan- 
tum computers might be able to solve problems that had no efficient solution on a classical 
or probabilistic Turing machine. He also presented an example oracle problerr0, now known 
as Deutsch's problem, that was solved more efficiently on a quantum computer than by 
any classical algorithm. In addition, four years later Deutsch [15] developed another, more 
popular, model for quantum computing known as the quantum circuit model. Nevertheless, 
neither Deutsch's problem nor its latter developed generalisation Deutsch- Jozsa's problem 
[16] proved to be practically useful, so the field continued to develop slowly. 

Peter Shor [57] shocked the computing world in 1994 when he presented two efficient 
quantum algorithms, for factoring integers and the 'discrete-log' problem, for which there 
was, and still is, no efficient classical counterpart. Researchers had been searching for an 
efficient factoring algorithm for over 20 years, and most people were confident that this al- 
gorithm did not exist as the most efficient classical algorithm discovered is the number field 
sieve [4T] that requires exponential time. In contrast, Shor's factoring algorithm requires 
only polynomial time. Shor's algorithm was also the first non-oracle quantum algorithm 
produced. Shor's algorithm was inspired by Simon's quadratic-time quantum algorithm 60 
that solves an oracle problem, which requires exponential time on a classical computei[f|. 
Nevertheless, all the focus was on Shor's factoring algorithm, because of its extremely im- 
portant applications in cryptography. All classical cryptography techniques, such as RSA, 
would be easily breakable with a quantum computer running Shor's algorithm. Therefore, 
it was this algorithm that spurred much interest and research in quantum computing. 

Even with this increase in research as a consequence of Shor's breakthrough, the only 
breakthrough since then has been Grover's search algorithm [29] developed in 1996. Grover's 
algorithm searches an n-element unstructured list in 0(^/n), compared to the classical 0(n). 
Thus, Grover's algorithm only provides a quadratic speed-up, but it is a major breakthrough 
due to its wide-spread application in search problems. Since 1996, the only new quantum 
algorithms have been variants of Shor's and Grover's algorithms. Possible reasons for the 
slow discovery of quantum algorithms are presented in Section [2.4.31 

2 An oracle problem is such that some properties of a function are trying to be discovered, whereby that function 
is given as a black box. The code of a black box function is unknown, but the function's results for given inputs 
can be calculated. 

3 Simon's algorithm, like Deutsch's, lacks practical application. 
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2.2 Preliminaries 



Depending on the reader's knowledge of linear algebra, computer science and quantum 
computing, some of the topics introduced in this section may already be familiar to the 
reader. Therefore, as these topics can be read independently, it is possible to just refer to 
those topics of personal interest. 

2.2.1 Dirac Notation 

Quantum states can be described as vectors, which are by convention expressed in the 
notation invented by Paul Dirac in 1958 |18| . The basic two elements of the Dirac notation 
are called bras and kets. Standard column vectors in a Hilbert Spacffl are represented by 
kets, such as \v) for a column vector v. The matching bra is a row vector, denoted (v\, 
which represents the conjugate transpose, or dual, of v. The conjugate transpose of v is 
defined as the transpose of the vector, in which each element is the complex conjugate of the 
corresponding element in v. These single vector representations can then be combined to 
represent such operations as the inner and outer product of vectors; for example, the inner 
and outer product of v and ip are written as (v\i/j) and |v)(y>| respectively. 

2.2.2 Tensors and Tensor Products 

In essence, a tensor is a geometrical identity that is the generalisation of n-dimensional 
vectors with their associated linear operators. Tensors can also be represented as multi- 
dimensional arrays. However, from a quantum computing standpoint, the tensor product is 
the most important operation that can be performed on tensors: the tensor product of v 
and ip is denoted v CS> if;. 

The tensor product is the most general bilinear operation: it is a generalisation of the 
matrix product operation, whereby all values contained in an operand are multiplied inde- 
pendently with all values in the other operand. Thus, the result of v ® ip has the following 
properties: 



The actual operation of the <S> operator is more clearly described by the following example: 



A complete treatment of the mathematics of tensor products is unnecessary here considering 
the scope of this paper; for a more complete treatment of tensor products refer to Hunger- 
ford's book [34] . The application of tensor products in quantum computing is discussed in 
Section 15X31 and 

2.2.3 Quantum Bit 

A quantum bit, commonly referred to as a qubit, is the basic unit of information in quantum 
computing: it is analogous to the bit in classical computation. Similar to a classical bit, 
a qubit has two computational basis states, sometimes referred to as eigenstates, usually 
represented as |0) and |1), that correspond to a bit's and 1 states respectively. However, 
unlike a classical bit, a qubit can be in a superposition of these two basis states. Due to this 
superpositional nature of a qubit it can be thought of as a 2-dimensional vector (in complex 
vector space) of length one with the above two basis states as orthonormal vectors: 



4 Chapter 16 of Hardy and Steeb's book I32| contains a good introduction to Hilbert Space. 



Rank(A B) 
Dimension(A <g> B) 



Rank(A) + Rank(B) 
Dimension(A) x Dimension(B) 
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Thus, a qubit in a superposition of basis states can be described as being a unit vector that 
lies between the two basis states. A qubit in an arbitrary state is therefore expressed as a 
linear combination of basis states: 

a|0)+/3|l) or 
where a, (3 are complex numbers and \a\ 2 + \f3\ 2 = 1 

Nevertheless, just like a classical bit, a quantum bit can be read to get its value, which can 
only be one of the basis states. However, the outcome is not deterministic as in classical 
computing, but rather probabilistic. Given the above expression for a qubit, the probability 
of the qubit being measured in each basis state is determined by the values of a and (3, which 
are referred to as amplitudes or probability amplitudes. The actual probability of the qubit 
being measured in a basis state is the square of the corresponding amplitude; for example, 
the probability of the qubit being measured as a 1 0} is \a\ 2 . Hence, the need for the condition 

M 2 + |/?| 2 = i. 

As the two amplitudes of a quibit in a basis state are 1 and (order irrevelant) it 
will have 100% probability of being measured in the state it is in and 0% probability of 
being measured in the alternative basis state; that is, a qubit in a basis state is measured 
dctcrministically. On the other hand, measurement of a qubit in a superposition of basis 
states changes the qubit into the basis state in which it is measured. The ramifications of 
this are that subsequent measurements of a qubit initially in superposition are deterministic 
and will yield the same outcome with 100% probability. Therefore, a qubit in a basis state 
can be thought of as a classical resource; in contrast, a qubit in a superposition of basis 
states can be thought of as a purely quantum resource |74j . 




2.2.4 Multiple Qubits 

Individual bits in classical computation combine through cartesian product, but quantum 
bits combine through tensor product [54] ■ Taking the simplest case of multiple qubits, two 
qubits, we can find the possible basis states as the tensor product of the two individual sets 
of basis states: 

{|0),|1)}®{|0),|1>} 
= {|0) (g> |0), |0) ® |1), |1) ® |0), |1) (g> |1)} or more concisely : 
{|00),|01),|10),|11)} 

The final concise version represents an extension of Dirac notation whereby \qoqi ■ ■ ■ q n ) 
represents the basis for qubits qo...q n . Thus, with k written as a binary number, a two- 
qubit quantum system can be described as: 




3 

where 8^ are complex numbers and S ^ \ 5kY 



Thus, an n-qubit quantum system can be written as X/fe-c) 1 ^1^)' where X/fc-o 1 l^ fe | 2 ~ ^ 
and \Sk\ 2 is the probability of the system being measured in the basis state | k). It is possible 
to measure a sub-set of the qubits in a multi-qubit system. It is also possible to determine 
such information as whether two qubits are equal without learning their value; however, these 
more complex measurements are equivalent to a transformation of the quantum system 
followed by a standard measurement to determine the basis state of each qubit [SI], and 
consequently it is common practise to also refer to standard measurements. 



2.2.5 Entangled States 

A suprising and non-intuitive aspect of quantum computing is that there are quantum states 
that can not be described in terms of their individual component qubits. These states are 
known as entangled states, and the individual qubits within entangled states are known as 
entangled particles. The reason for the use of the term 'entangled' should become obvious 
in the next paragraph. 
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We established in the previous section that single qubits combine using the tensor product 
operator. Therefore, two qubits in superposition combine as follows: 

I gi } ® |q 2 > = (ai|0)+ ( 3 1 |l))®(Q 2 |0>+/3 2 |l» 

= aia 2 |00) +ai/3 2 |01) +a 2 /3i|10) +ft/3 2 |H> 

However, there are entangled states that can not be described by the above formula. The 
canonical examples are the Bell states: 

hM = ^^ hM = ^P 

In order to write Voo in terms of the above two-qubit equation ai/3 2 must equal 0. This 
subsequently implies that either cvitv 2 = or fiifh = 0; however, it is clear that neither of 
these terms are zero. Thus, ipoo is an entangled state (the other Bell states are shown to be 
entangled with similar reasoning). An example of how a non-entangled state can be expressed 
in terms of its component qubits is the state s ^ rt2 (1 10) + 1 11)) = 1 1) <8> s ^ rt2 (1 0) + 1 1)) 

Another way of viewing entangled states is that measurement of one qubit has an effect 
on the other qubits. Take ipoo as an example again, if the first qubit is measured as |0) 
then we know that the second qubit will also be measured as j 0). In contrast, measurement 
of the first and second qubit in our non-entangled state will always yield |0) for the first 
and a fifty-fifty chance of either basis state for the second, regardless of the measurement 
of the measurement order. Some famous research explaining this observation includes that 
of Einstein, Podolsky and Rosen: an accurate and concise summary of this work is given 
by Rieffel and Polak [54]. It is clear that there is no classical counterpart of entangled 
states. However in quantum computing, it is important to realise that entangled states 
are treated no differently from any other quantum states. Nevertheless, they are extremely 
important and some very interesting applications of entangled state have been found such as 
dense coding, quantum teleportation (see Section ^. 7. l|l and secure quantum key distribution, 
which are discussed by Rieffel and Polak [SI]. A detailed analysis of entangled states, their 
importance and a C++ implementation, is provided by Hardy and Steeb [31] . 

2.3 Quantum Computers 

Quantum computers still remain abstract machines despite the large amount of money and 
time that has been invested into building them. Quantum computers with a few qubits have 
been developed; however, a scalable general- purpose quantum computer is yet to be built. 
The largest quantum computer developed thus far consists of seven qubits [69], which was 
built by a group of scientists from Stanford University. It was constructed using Nuclear 
Magnetic Resonance (NMR) and tested using Shor's factoring algorithm J7SJ. In February 
2007, the Canadian company D-Wave publicly demonstrated 'Orion', which they claim is 
a 16-qubit quantum computer. However, details of Orion's inner workings are unknown 
as no academic papers have been published about OriorQ Thus, there is academic spec- 
ulation about whether Orion is truly a quantum computer. There have been many other 
approaches to building a quantum computer, including optical photon spins, quantum dots, 
cavity quantum electrodynamics and ion traps. NMR and ion traps are the most advanced 
approaches to date; however, all of these techniques have had scaling problems. Therefore, it 
is commonly argued whether or not quantum computers of a significant scale (greater than 
100 qubits) will ever be practically realised. Pellizzari [48 presents a general overview of the 
requirement for building a quantum computer, including its challenges, various approaches 
and future outlook; however, this article was published in 1998 and consequently does not 
include the new NMR approach. Hardy and Steeb [35] provide a concise summary of NMR, 
and Nielsen and Chuang [47] provide a more detailed look at all of the above mentioned 
approaches. 

There are certain characteristics that an operational quantum computer must have re- 
gardless of the approach adopted for its construction. Different publications outline various 
characteristics that are somewhat similar 32 , 42, 47 . These characteristics have been com- 
bined and summarised below: 

5 One of many news articles about Orion can be viewed at 
http: / / arstechnica.com/articles/paedia/hardware/quantum.ars 
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Qubit Representation There must be a representation of a qubit (dynamic two-state 
object) such as the up and down spins of a proton. 

State Preparation The system must reliably start in a known initial state. It is sufficient 
and common to start with all qubits in the same state (conventionally |0)); thus, this 
requirement has not been a major issue. 

State Transformation There must be a mechanism to efficiently transform a quantum 
system into other states following the rules of quantum mechanics. The most common 
way this has been implemented is by the use of quantum gates, which will be discussed 
in Section EH] 

Measurement There must be an efficient and reliable way to measure qubits in the system. 

System Isolation The system must be isolated from the environment to prevent the su- 
perposition of states from decaying quickly in the common environment, which is a 
phenomenon known as decoherence or quantum noise. This is a problem as algorithms 
may not be able to be applied in time if the states decay too quickly. This issue of 
decoherence is a major hurdle for the actualisation of a quantum computer 54, 37, 59 . 
Zurek 84 provides a comprehensive and simple discussion of decoherence. 
It is almost inevitable that complete isolation will not occur and there will be some 
decoherence. As a partial solution, research is being undertaken on quantum error- 
correcting codes (QECCs). Since the first 9-qubit QECC was presented by Shor in 1995 
[58] research and development in this area has been making noticeable progress; the 
most recent work presents the most time-efficient QECC (although not the most space- 
efficient) [23]. Despite progress being made in developing quantum error-correcting 
code, this method will always come at a cost of extra qubits, which consequently 
increases the number of qubits needed to implement any given algorithm. This is the 
reason QECC is solution for only a small amount of decoherence. Nielsen and Chuang 
[47] present a comprehensive section on the area of quantum error-correction, which 
has been summarised by Landry [37]. Furthermore, Hardy and Steeb [35] discuss the 
latest algorithms for developing quantum error-correcting code. 

Whether a scalable quantum computer will ever be built obviously has huge impacts 
for the field of quantum computing, as well as classical computing. Nevertheless, quantum 
computing is currently a popular research field and this paper simply views a quantum 
computer as an abstract machine. Furthermore, while this paper does not discuss quantum 
mechanics directly, most of the basic postulates of quantum mechanics will be discussed 
within Section [5] of this paper. For a structured coverage of the exact postulates of quantum 
mechanics refer to Nielsen and Chuang's book [47] . 

2.4 Power of Quantum Computing 

After more than 50 years of using the classical physics paradigm (which the classical Turing 
machine embodies) to build a theory of (standard) computation, quantum physics provides a 
different paradigm that is arguably more powerful than standard computation [4]. The pos- 
sibility of harnessing the power of quantum parallelism (discussed in the following section) 
and the identification of entangled states that are without a classical counterpart (discussed 
in Section I2.2.5P were the first indications that quantum computing could allow faster in- 
formation processing than classical computing. Furthermore in 1985, Deutsch proposed 
that his quantum Turing machine could not only simulate quantum systems better than 
classical methods, but it may also be able to solve classical problems significantly faster 
than classical, and possibly probabilistic, Turing machines. This essentially suggested a vio- 
lation in the strong Church- Turing thesis that any algorithmic process can be simulated on 
a classical Turing machine (CTM) in at worst a polynomial slowdown. 

2.4.1 Quantum Parallelism 

Quantum parallelism is the term used to describe the potential parallel computing power 
of quantum computing. It is intuitive that an increase in time efficiency can be gained 
by using parallel processors. Furthermore, in classical computation, an exponential time 
efficiency increase requires an exponential increase in the number of processors or physical 
space. However, as qubits can represent a superposition of two different states, in quantum 
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computation, a linear increase in physical space generates an exponential increase in paral- 
lelism, and consequently an exponential increase in time efficiency. This is what is known 
as quantum parallelism [16] . As shown in Section 12.2.31 a qubit can be represented as a 
two-dimensional vector. In classical physics, n two-dimensional vector objects form a 2n- 
dimension vector space. However due to quantum parallelism this is not the case in quantum 
computing in which n qubits form a 2"-dimension vector space. This exponential growth 
implies a possible exponential increase in the information processing speed of quantum com- 
puters over classical computers. The concept of how quantum parallelism is applied to create 
parallel computation is conveniently thought of as shown below in a two-qubit example: 

(!) • - 

All Inputs <g> One Transformation 

A set of input qubits in a superposition of basis states essentially has all the possible inputs 
encoded in them. Therefore, a single transformation can be applied to generate a set of 
output qubits in a superposition of basis states, which represents all possible outputs. Thus, 
all output states have been computed and assigned a probability of measurement on a 
quantum computer in the same time it takes to compute the output for a single input state 
on a classical computer. 

Although all outputs can be produced using only one transformation, we know from 
Sections 12. 2. 31 and 12. 2. 41 that only one output can be extracted upon measurement. Further- 
more, the particular output that will be extracted is unknown as measurement is probabilis- 
tic. Thus, a quantum algorithm must manipulate quantum parallelism so that the desired 
result is extracted when measured. This is a difficult and non-intuitive task that no classical 
programming technique can solve. The two best-known approaches are [54] 

• Measure common properties of all the outputs. This approach is used in Shor's factor- 
ing algorithm whereby the period of the outputs is measured 

• Increase the amplitudes of the basis states of interest to increase the probability that 
they will be measured. This technique is used in Grover's search algorithm. 

In addition, there are other restrictions on quantum computing that are discussed in the 
next section. 

2.4.2 Unitary Restriction 

With the exception of measurement, all transformations in a quantum system must be uni- 
tary. This is specified in Schrodinger's equation, which governs the dynamics of a quantum 
system which is not being measured. To be precise, Schrodinger's equation enforces that 
quantum systems must change state in a way that preserves orthogonalitjlf], and the only 
transformations that maintain orthogonality in a complex vector space, such as quantum 
systems, are unitary transformations. The transformation matrix M is a unitary matrix if 
MM* = I, where M* is the conjugate transpose of M. A unitary transformation can be 
thought of as a rotation in complex vector space [54] : this makes sense as vector rotations 
will maintain the angles between all vectors and preserve orthogonality. 

An important consequence of quantum transformations being unitary is that they must 
be reversible, as unitary matrices are invertible. For computing to be reversible, the current 
state must uniquely determine the previous state, for all computational states 17" . However, 
this is not a problem as Bennett [7] showed that all classical computation can be transformed 
into reversible computation without a polynomial expenditure in time or space. Another im- 
portant consequence of unitary transformations is that their linearity implies that quantum 
states can not be cloned, which is known as the no-cloning theorem. An easily understand- 
able version of the original proof of this theorem by Wootters and Zurek [80] is given by 
Rieffel and Polak [SI]. It is important however to note that this theorem only applies to 
qubits in an unknown state: qubits in a basis state or a known state of superpositiorQ can 
be cloned, but qubits in an unknown state can not be cloned as the information about their 
amplitudes will be lost. 

6 To be orthogonal means to be linearly independent: a precise treatment of the orthogonality condition is 
given by MathWorld [available at |http:// mathworld, wolfram.com/OrthogonalityConditi on. html] . 
7 A known state of superposition is where all the amplitudes are known. 
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2.4.3 The Potential 



There is a well-established view that quantum computing will only yield an exponential 
speedup in problems whose structure avoids the need to process exponentially many cases 
[12] . That is, a brute force approach to NP-complete problems will not succeed with the 
aid of quantum parallelism; Fortnow and Rogers [24[ also firmly question whether quan- 
tum computing is more powerful than classical computing, but they still maintain it is a 
potentially powerful model of computation. However, whether quantum computation can 
efficiently solve NP-complete problems in polynomial time remains an open question [54] . 
With the discovery of a few quantum algorithms there is no doubt that certain problems can 
be solved more efficiently on a quantum computer than on a classical computer. The number 
and diversity of problems that can be solved more efficiently on a quantum computer is still 
unknown and the subject of much current research. 

In reality, it is not extremely difficult to develop a working quantum algorithm. How- 
ever, quantum algorithms are only of interest if they are more efficient than their classical 
counterparts. This is one of the main reasons that more quantum algorithms have not been 
discovered. It is simple enough to simulate a classical algorithm on a quantum computer; 
however, producing a quantum algorithm that is superior to its classical counterpart requires 
the use of truly quantum effect^j, which is extremely complex. Pessimists speculate that 
the lack of discovery of quantum algorithms is due to a lack of application of quantum com- 
puting. However, it seems more likely the reasons are the non-intuitive nature of creating 
quantum algorithms to the classically trained world and the fact that only superior quantum 
algorithms are of interest. Nevertheless, it has not been conclusively proven that quantum 
computation is more powerful than classical computation (this issue will be discussed more 
in the following section). The potential of quantum computation is however extremely high 
as there may even be an exponential speed gap between a quantum and classical computer 
[27] . as is the case for Shor's factoring algorithm [57] . 

2.5 Complexity of Quantum Algorithms 

Quantum computing is a unique type of probabilistic computing, which has caused the 
creation of new complexity classes [54]. The most interesting quantum complexity class 
is bounded-error quantum polynomial time (BQP), which is the quantum equivalent of 
the classical bounded-error probabilistic polynomial time (BPP) complexity class. BPP 
contains all languages for which there is a probabilistic Turing machine such that it gives 
the right answer with bounded probability, which means greater than or equal to | of the 
time. Note that the probability | is simply a norm as it can be replaced with any number 
between \ and 1 without altering the complexity class: if the probability is set higher than 
i any previously BPP algorithm can simply be rerun to amplify the probability [SS]. Using 
this definition of BPP and substituting 'quantum Turing machine' for 'probabilistic Turing 
machine' we obtain the definition of the complexity class BQP. 

It is a significant result in quantum complexity theory that Bernstein and Vazirani proved 
that BPP C BQP C P SPACE [5], but the open question of extreme importance in 
quantum computing is whether BQP is strictly larger than BPP. It is important to note 
that proving BPP C BQP implies that BPP C PSPACE, which is not currently known. 
Thus, proving that quantum computers are more powerful than classical computers would 
also represent a breakthrough in classical computing complexity theory. However, it also 
indicates that proving BPP C BQP is probably extremely difficult. Simon's quadratic-time 
quantum algorithm [BO] to solve a problem with known classical complexity of exponential- 
time suggests that BQP is strictly larger than BPP, but it is not a rigorous proof as it 
only shows BPP C BQP for oracle problems [59[ . Grover's search algorithm [25] also puts 
forward a case that BQP is strictly larger than BPP; however, it remains an open question 
[39] . That is, it has not yet been proven that quantum computers have more capabilities 
than classical computers. 

Further details on quantum complexity classes are outside the scope of this paper, but 
more information on the BQP complexity class is presented by Fortnow and Rogers [24] . 
and [5J [77] detail quantum computing complexity in time and space respectively. 

8 The term quantum effects is used as an umbrella term to encompass effects such as superposition and entan- 
glement that exist in quantum computing, but not classical computing. 
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2.6 Models of Quantum Computation 



The first quantum computation model was the quantum Turing machine (QTM) introduced 
by Benioff 6 and Deutsch [14], which can efficiently simulate every classical Turing machine 
(CTM) with a polynomial-bounded overhead ¥l\ . The QTM is also the commonly used basis 
for defining universality. In essence, the QTM has three deterministic CTM tapes where the 
two extra tapes are needed to make the computation reversible (the set of tape states also 
changes to reflect the nature of quantum computing). The QTM is described in more detail 
by Gramfi [27]. 

There are also many other mathematically equivalent models of quantum computation, 
just as there are many different models of classical computation. Moreover, there are other 
potential models of quantum computation that have not been explored [59] . A quantum 
state machine [30] has been proposed that is an extension of classical finite state machines 
in which amplitudes are added to transitions to represent quantum parallelism; a universal 
quantum cellular automata model [73] is another quantum computing model. Recently, a 
universal quantum lambda calculus [74] has also been suggested that is based upon the 
classical lambda calculus. Essentially, a reversible lamdba calculus is established with the 
use of minor reduction steps, and the quantum parallelism is captured as subexpressions 
representing superposition that can not be reduced/observed (no-cloning theorem). 

This section of this paper focuses on the quantum circuit model, due to its popularity 
and simplicity. The quantum circuit model, also known as the quantum gate array, was 
introduced by Deutsch in 1989 [15] . The circuit model is mathematically equivalent to the 
QTM 83 , and is the most popular model of quantum computation. The reason for the 
popularity difference is that the QTM suffers from the same complexities as the CTM; it 
does not satisfy a simple algebra and it can be cumbersome to use as it requires 'word at a 
time' thinking while keeping track of control variables such as tape states [71]. In contrast, 
quantum circuits are easier to understand as they are graphical diagrams, which can also be 
manipulated algebraically. Furthermore, other models such as the quantum lambda calculus 
are too new to be widely popular. The quantum circuit model will be explained in more 
detail in the following sub-section. 

2.6.1 Quantum Circuit Model and Quantum Gates 

Similar to classical circuits, the quantum circuit model consists of wires, and gates that act 
upon wires. In the quantum circuit model, each wire represents a qubi10 and each individual 
gate has the same number of input qubits as it has output qubits (due to the reversibility 
of quantum computing). Thus, an advantage of this model is that the diagram itself can 
be essentially thought of classically, with the data on the wires representing the quantum 
nature of the computing. As the number of qubits (wires) is required to be constant, an 
important restriction of this model is that a quantum circuit can only compute a function 
with a specific length domain. Therefore, for functions with arbitrary length inputs a family 
of quantum circuits is required; that is, a quantum circuit for each input length is required. 

The complexity of a quantum circuit depends on the total number of gates and qubits 
(wires). However, as gates are commonly restricted to act upon one, two or three qubitsH 
the number of gates in a quantum circuit is both a reasonable and the usual measure of 
complexity. This restriction does not affect the universality of quantum circuits; for example, 
the set of all one-qubit gates and the CNOT gate (explained later in this section) are universal 
for quantum computing. Further note that an interesting three-qubit gate is the Toffoli gate 
that is in fact universal for quantum computing as presented by Deutsch [15] . In addition, 
it has more recently been shown that a large number of two-bit quantum gates are universal 
[3J- An excellent review of quantum gates and the respective universality issue is contained 

in @. 

Quantum gates are in fact unitary transformation matrices (see Section I2.4.2[l , and a 
matrix that acts upon n-qubits will be a 2" x 2™ square matrix due to quantum parallelism. 
Quantum parallelism (see Section \2. 4. l|l is the reason for a linear increase in the number of 
qubits causing an exponential increase in the size of a quantum system. To familiarise the 

9 By convention a qubit begins in a basis state, usually assumed to be |0). 

10 There is some inconsistency in the literature about the maximum number of qubits that a gate should act 
upon; for example, Shor [59] states the maximum should be two qubits while Spector et al. |63] state that the 
maximum is a few qubits. 
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reader with both the gates and the quantum circuit modef^], some important elementary 
gates are detailed below using the quantum circuit model. 

Three important single-qubit gates are the identity (I), negation (X), and Hadamard (H) 
transformations, which are detailed belotf^l: 



X 



H 



1 

1 

1 

1 



>2 



1 
-1 



10} 



10} 



10} 



|0);|1) 



|i>;|i> 



|i} 



10} 



V2 



(|0) + |1));|1) 



1 



|0>-|1» 



It is important to notice that the Hadamard gate creates an equal superposition of states, 
and when applied to n-qubits it creates a superposition of all 2" states (this is the reason 
for its prevalence in quantum algorithms) . The workings of these gates can be easily verified 
with linear algebra; for example, the negation transformation on |0): 



1 

1 



1x0+0x1 
1x1+0x0 



These single-qubit transformations are represented in the quantum circuit model by an 
appropriately labeled box as shown below: 





X 











H 









The controlled-not (CNOT) is an extremely important two-qubit gate. Its importance 
comes from the fact that the results of a CNOT transformation can not be expressed as the 
tensor product of two qubits; thus, a CNOT gate can be used to create entangled states. 
The workings of the CNOT gate are relatively simple: it swaps the target qubit (represented 
by a cross) if the control qubit (represented by an open circle) is |1). The mapping and 
quantum circuit notation for the CNOT gate are as shown. 

1 00) -> 1 00) 1 01) -» 1 01) 
1 10) -> |11) |11) -> 1 10) 



The universal 3-qubit Toffoli gate is essentially a controlled-controlled-not gate, where 
the target qubit is only flipped if the two control qubits are |1); It is represented in the 
circuit model as follows. 



-e- 



-o- 



11 Variations in the quantum circuit notation of different gates exist, but they are only minor and should be 
understood by a reader who has read this paper. 

12 Due to linearity, transformations are fully specified by their effect on the basis states |54| . 
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2.7 Example Problems and Algorithms 



The algorithms that will be presented in this section are quantum teleportation, and an 
overview of Grover's algorithm. These algorithms have been chosen for their ease of under- 
standing and their ability to convey important quantum algorithm concepts. Pittenger's 
book [3?5] contains excellent, succinct treatments of the following quantum algorithms: 
Deutsch-Joza's, Simon's, Grover's, Shor's and the Quantum Fourier Transformation. 

This section will also include definitions of Deutsch's and the scaling majority-on prob- 
lems. Quantum algorithms to solve these problems have been evolved by Spector et al. using 
genetic programming (see Section [4. 3. 3|l . 



2.7.1 Quantum Teleportation Algorithm 

Quantum teleportatiorF^l. uses two classical bits to transfer one quantum bit of information 
(|gl)) from A(lice) to B(ob) without ever being anywhere in betweerP^I. This does not 
violate the no-cloning theorem as the unknown qubit is teleported not copied; that is, after 
the quantum teleportation the qubit no longer exists with Alice. The quantum circuit model 
for teleportation is as shown below in Figure [1] 



A 
B 



-e- 



ii 



-e- 



ii 



-X- 



Figure 1: Quantum Teleportation Algorithm 

The first step creates Bell's first entangled state IV'oo) = 00 ^j 11 ' 1 (discussed in Sec- 
tion [223} by putting the first qubit into equal superposition using the Hadamard gate and 
then modifying the second according to the value of the first using the CNOT gate. Now, 
given that |<jl) is unknown, it can be represented as a\0) + f3\l), and the quantum state 
after the first two transformations is 

-J=(a|0)(|00) + |ll))+/3|l)(|00) + |ll))) 

Alice and Bob are each given one of these two entangled particles and the final two trans- 
formations are performed as shown, and the next CNOT is applied to make the quantum 
state 

-^(a|0)(|00) + |ll))+/3|l)(|10> + |01))) 

Thus, the information about \ql) is now contained in the entangled pair, and the final 
Hadamard transformation is performed to yield 

|(a(|0) + |1» (]00) + |11)) +/3(|0) - |1»(|10) + |01») 
= i[|00) (a|0)+/3|l» + [01) (a|l)+/3|0)) + |10) (a|0) - /3\ 1)) + 1 11) (a\l) - p\0))} 



Alice then measures the two qubits she has and sends the result encoded in two classical 
bits {00,01,10,11} to Bob. Using the previous equation, Bob knows the state of his qubit 
in terms of |ql). Therefore, he can apply a simple one-qubit gate (if needed) to convert his 
qubit to \ ql) as defined below (using the previous equation): 

Alice's Result Bob's Qubit Transformation Required 
1 00) a\O)+0\l) none 
1 01 > a\l)+P\0) negation(X) 



10) a\0)-/3\l) 



1 

-1 



13 This explanation of the quantum teleportation algorithm is based on the explanation given by Landry [37] ■ 
14 Quantum teleportation of one qubit has been realised experimentally pill . 
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|H) «|1>-/3|0) ( I _ j 
For example, in the case of |01), after transformation Bob's qubit equals 

*(«ii>+0o» = (; ;)(; 

= a|0)+/?|l) 
= l«l> 

Thus, an unknown qubit in a superposition of states can be transported from A to B with 
only two classical bits of information. 



2.7.2 Grover's Algorithm 

Grover's algorithm searches an unstructured n-element list in 0(y/n) time compared to its 
classical counterpart which is 0(n). An overview of the algorithm in the quantum circuit 
model is shown below in Figure [2] (based on Figure 6.1 in Nielsen and Chuang's book 1471 ). 



qubits 



10] 



oracle 
workspace 









G 




G 














































Figure 2: Grover's Algorithm 



The gate labeled H is a shorthand notation for applying a Hadamard gate to the n- 
qubits. The gate labeled G is known as the Grover iteration or Grover operator. Essentially, 
the Grover operator increases the amplitude of the basis states that are being searched for 
and decreases the amplitudes of the other states. This is done by rotating the current state 
vector (of the n-qubits) towards the superposition of all solutions to the search problem 
(|/3)) as shown in Figure [3] The Grover operator consists of 

1. Applying the oracle search function (using the oracle workspace) to the n-qubits. 

2. Applying Hadamard transformations to all n-qubits 

3. Performing a conditional phase shift 

4. Applying Hadamard transformations to all n-qubits 

The initial oracle application reflects the current state of the n-qubits about \a) to move 
it from its initial state 1 to state 2, as shown in Figure [3] The remaining three operations 
(Hadamard, conditional phase shift, Hadamard) then reflect the state 2 vector about state 1. 
Overall, this achieves a rotation of 9 towards \ j3) (state 3). Repeating the Grover operator 
0(i/n) times rotates the state vector very close to |/3), the superposition of all solutions. 
Hence, when the state is measured with reference to the computational basis it reveals 
a solution to the search problem with extremely high probabilitj03- This is the essence 
of Grover's algorithm; an excellent detailed analysis of Grover's algorithm and its recent 
variations is given in chapter 6 of Nielsen and Chuang's book |47j . 



2.7.3 Deutch's and Scaling Majority-On Problems 

Deutsch's problem is determining whether a given oracle, or black box, function is uniform 
or balanced, given that the oracle must be either one of these. The uniform property 

The exact probability of observing a solution depends upon the number of solutions in the search space. 
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Figure 3: The result from applying the Grover operator: modified from Figure 6.3 in Nielsen 
and Chuang's book [J7]. The three numbers in boxes indicate the initial (1), intermediate (2), 
and final (3) state vectors when the Grover operator is applied. \(3) is the superposition of all 
solutions to the search problem, and |a) is a state orthogonal to \j3). All vectors in this diagram 
are unit vectors, but \a) and \/3) are shown to be longer to improve clarity. 

requires an oracle to always return or always return 1; the balanced property requires 
an oracle to return an equal number of Os and Is over all possible inputs. The scalable 
majority-on problem is an extension of Deutsch's problem where the oracle is an arbitrary 
boolean function and the problem is to determine whether the majority of the outputs are 
f . Quantum algorithms to solve these problems have been evolved by Spector et al. using 
genetic programming (see Section H. 3. 3 p - 

2.8 Further Reading 

This section presents the major summary references in the field of quantum computing that 
complement the specific references given in previous subsections. However, most of the 
references given here have been previously referenced in this paper. 

There are currently a large, and ever increasing, number of publications on quantum 
computing, and consequently there is a lot of duplication and work of varying standards. 
For issues on quantum computing, Nielsen and Chuang's book 47 appears to be the most 
comprehensive and well-structured publication, which is considered to be the most significant 
reference. A very good introductory reference is provided by Landry [37], who summarised 
the major areas in Nielsen and Chuang's book. Another excellent and recommended review 
reference is Rieffel and Polak's work |54J . which contains a more comprehensive and detailed 
review of quantum computing than Landry's paper [37] . Furthermore, in contrast to [371154] . 
Gramfi [27] provides a summary of quantum computing from a quantum Turing machine 
perspective, rather than the quantum circuit model. However, many (although not all) of 
the topics included in [27\ I37[ 154] have been covered in this paper. 

There are two excellent quantum computing books that also contain detailed informa- 
tion on how to actually simulate quantum algorithms on a classical computer: Williams 
and Clearwater's book [79] comes with MathematicaF^l notebooks that simulate well-known 
quantum algorithms such as Shor's factoring algorithm, and Hardy and Steeb's book [32] 
contains Java and CH — h code for some simple quantum simulations such as generating en- 
tangled qubits. 

The majority of the articles referenced in this quantum computing section are freely 
available at the Los Alamos preprint server: http://xxx.lanl.gov/archive/quant-ph This 
site also provides an excellent place to search for old and new articles within the field of 
quantum computing. 

16 Mathematica is a comprehensive mathematical software package, details are available on the their website: 
http: / /www. wolfram.com/products/mathematica/index. html 
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3 Introduction to Evolutionary Algorithms 



Optimisation problems can be characterised by two sets of parameters: feedback param- 
eters to optimise according to a target solution, and free parameters to modify in order 
to approach the desired solution [1U| . The optimisation algorithm alters the free param- 
eters while controlling the quality of the solution by the feedback parameters; different 
optimisation techniques perform this search in different ways. When traditional optimisa- 
tion techniques are used to search vast, complex and unknown spaces there are extremely 
complex constraints and multimodal problem£3 [32]. Thus, traditional optimisation tech- 
niques are not well-suited for these types of problems, so alternative approaches have been 
researched. Evolutionary Algorithms (EAs) are one of the alternative methods that have 
gained significant popularity as general-purpose optimization and search tools [32]. EAs are 
probabilistic search algorithms that are heavily based upon Darwinian evolution, as a proxy 
for the process of species evolution in nature. The central concept of Darwinian evolution 
is that individuals in a population have heritable characteristics that influence their prob- 
ability of producing offspring, that is, future generations. EAs have extended this theory 
slightly to state that characteristics of 'better' individuals will increase their likelihood to 
produce offspring, which is a variation on Darwin's 'survival of the fittest' principle [13] . 
The idea behind this is an attempt to converge to the 'best' individual, which essentially is 
the paradigm of search and optimisation. 

EAs were initially used as optimisation tools for engineering problems, and were devel- 
oped independently by several computer science researchers in the 1950s and 1960s 25 . 
Since then, the number of applications for EAs has become diverse and has grown at a fast 
rate [10] : for example, there are EA applications in financial forecasting, predicting protein 
structure, predicting the primeness of numbers and in developing computer programs. In 
contrast to the application aspect of EAs, despite much theoretical research there has been 
modest progress in EA theory over the last 20 years compared with that of neural networks, 
another biologically motivated form of computation [5]. Nevertheless, various types of EAs 
have been developed and the types presented in the following subsection now form the back- 
bone of the EA field [46]. Regardless of what type of EA is used the basic elements are 
almost identical |2]; these common elements are [39] 

• Populations of individuals representing solutions to the problem at-hand, which allows 
parallel searching. 

• Ways to manipulate solutions, which can be either 

— Mutation (inspired by the biological process of the same name) operators, which 
implement innovative change 

— Recombination/Crossover (inspired by the biological process of gamete production 
and sexual reproduction) operators, which implement conservation of characteris- 
tics through rearrangement 

• A measure for determining the quality of a solution, usually referred to as a fitness 
function 

• A method of selection that uses the fitness function to select individuals for the next 
generation 



3.1 Types of EAs 



Evolutionary Algorithms 



Evolutionary Programming 



Genetic Algorithms 



Evolution Strategies 



Genetic Programming 



Figure 4: Relationships between different types of Evolutionary Algorithms 



17 Multimodal problems are problems that arise in cases with a large number of locally optimum solutions. 



14 



EAs and Genetic Algorithms (GAs) are terms that are sometimes incorrectly used inter- 
changeably. However, EA is an umbrella term that includes all algorithms that incorporate 
the idea of Darwinian evolution, and GAs are just one type of EA. The different types of EAs 
are shown in Figure [3] GAs, invented by Holland in the 1960s [33j, are the most prominent 
EAs 32 , with the other major types of EAs being Evolutionary Programming (EP) and 
Evolution Strategies (ES). Essentially, EP and ES operate on and change the phenotype 
(observable properties) of individuals, while GAs operate on the genotype (genetic construc- 
tion) of individuals [22]. To further explain the concept of phenotype, a change in phenotype 
can be viewed as a change in the behavior, physiology or morphology without altering the 
genetics. In addition, a further difference between ES and GAs is that ES places an em- 
phasis on mutation [10] . whereas GAs place a higher weighting on recombination. However, 
having stated the differences between these types of EA, it is important to note that these 
distinctions are not strict, as overlap does occur [TU]. For example, GAs have been applied 
to real numbered genomes [70], which overlaps into the ES field. For a more thorough review 
of EP and ES, which also includes a discussion on GAs, please refer to Back's book [T]. 

In 1992, Koza introduced a new type of EA known as Genetic Programming (GP) 36 ; 
Bornholdt 10 suggests that the reason GP arrived late into the EA field may have been due 
to the need for greater computing power. GP is a type of GA whereby the search space is 
reduced to solely include computer programs. GP techniques have been valuable in evolving 
structures other than computer programs, but the fact that individuals in a GP model are 
computer programs is the most defining feature of GP [61]. Nevertheless, GP has extended 
the idea of genotype manipulation from GAs to include variable length chromosomes; that 
is, the representation of individuals in a population can be of varying size. 

3.2 GA and GP Algorithm Structure 



generate 

initial 

population 



t 

stait 



evaluate 
fitness function 



new 

population 
created 



Are termination 
criteria met ? 




Darwinian Selection 



Genetic Modification 
(Recombination / Mutation) 



optimised 
individuals 



Figure 5: A general flow diagram of GAs and GP algorithms, modified from Figure 3.1 in Leier's 
thesis P] . 

The overall structure of GAs and GP algorithms (GPAs) is shown in Figure [5] Initially, a 
starting population is generated; this is usually done randomly so that it results in a diverse 
population that simulates all naturally occurring populations. However, before an initial 
population can be created there must be a way to represent individuals. Traditionally, as 
supported by Holland [33], individuals, or solutions, were encoded as binary strings; modern 
GAs (and consequently GPAs) use tailored encoding to suit the problem being solved |45j . 
The number of individuals in the initial population is another parameter that must be set 
according to the specific problem at hand. 

After an initial population is created the evolutionary loop begins. The evolutionary 
loop consists of the following repeating sequence: 

1. Evaluate the fitness of each individual in the population. The fitness, or quality, of 
each individual (x) is measured by a predetermined fitness function (f(x)) with a real 
number codomain, that is, it returns a real number. 
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2. Select individuals based on fitness levels to be parents of the next generation, as in 
Darwinian Selection. The three main selection methods used in GAs and GPAs are 

(a) Fitness-Proportional selection, where individuals are selected using a roulette wheel 
concept in which each individual (i) has a sector of size proportional to its fitness 
(/;). Assuming that a lower fitness function corresponds to a superior fitness, this 
means each individual has a probability of being selected equal to 

Ei(i-/«) 

(b) Ranking selection, which simply involves selecting the best x (predefined param- 
eter) number of individuals, according to a predefined function based on their 
fitness. 

(c) Tournament selection, which involves randomly selecting a number (predefined 
parameter) of individuals to compete in the tournament. The losers of the tour- 
nament, based on fitness, are then excluded from selection, and all remaining 
individuals are selected. This is however, a much weaker form of selection than 
the ranking selection. 

3. Generate the next population through genetic modification. The mutation operator is 
applied to a single individual which results in a random change to its representation 
(genome) at one or more positions; it is applied with a certain low probability, set 
as a parameter. The crossover operator is applied to two individuals (parents) and 
produces a new individual (child) for the next generation, which is a combination of 
the representations (genomes) of the original individuals. The crossover operator is 
also applied with a probability set as a parameter, but this probability is usually much 
higher that the mutation probability. Thus, cross-over can obviously be performed 
in many ways and depends on the representations of the individuals. Two possible 
approaches are (1) using the beginning of the representation from one parent and the 
end from the other parent, or (2) reusing one parent's representation with a middle 
section exchanged for a section of the other parent. One common classification of 
different crossover operators can be separated into the following two categories: 

• Fixed-length crossover, where the children have the same length representations 
as their parents. Assuming that individuals are a sequence of instructions, an 
example of fixed-length crossover is taking the first 4 instructions from parent 
A (length 6) and the last 2 instructions from parent B (length 6), which would 
produce a child that had 6 instructions, the same as its parents. 

• Variable-length crossover, which simply means that the result of the crossover 
operation may result in a different length child. For example, taking the first 4 
instructions from parent A (length 6) and the last 5 instructions from parent B 
(length 6) produces a child that had 9 rather than 6 instructions. 

This loop terminates when sufficiently optimised individual£f| have been evolved or other 
predetermined termination criteria are met, such as the maximum number of generations 
have been evolved. 

3.2.1 Further Reading 

For futher information on GAs, Bornholdt 10 presents a succinct review, and Hardy and 
Steeb [32] provide a good review of the implementation of GAs, which includes many Java 
and C++ code examples. In addition, a very comprehensive analysis of GAs from a mathe- 
matical standpoint is given by Vose [76]. More detail on GP in general is covered in Koza's 
books [36] [35]. 

4 Evolving Quantum Algorithms 

Genetic Programming (GP) is defined as developing algorithms based on Darwinian evolu- 
tion to automatically generate computer programs. Thus, as quantum computing is a formal 
computational model, GP can be used to evolve quantum computing programs (represented 

The desired level of optimisation is set as a parameter in terms of the fitness function. 
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as quantum algorithms); this field will be hereafter referred to as 'Evolving Quantum Algo- 
rithms'. Unlike other applications of GP, there has been relatively little work done in the 
field of evolving quantum algorithms [55] . 

As was stated in Section 12.4,31 it seems that the major reasons for a lack of discovery 
of quantum algorithms are that they are difficult and non-intuitive to generate, and that 
only faster than classical quantum algorithms are of interest. Therefore, it would seem 
appropriate that we look to computers to search for quantum algorithms. However, as 
the search space of quantum algorithms is vast, complex and reasonably uncharted, using 
traditional search techniques is problematic (see Section [3|. Hence the motivation to use 
GP to evolve quantum algorithms, with the hope that GP's power to randomly search vast, 
complex and unknown spaces can discover many superior quantum algorithms. 

4.1 Simulation Limitations 

In order to evolve quantum algorithms the fitness of different individuals (quantum algo- 
rithms) must be assessed; the fitness of a solution is assessed by running the algorithm. 
However, as quantum computers have not yet been built (Section 12. 3p . quantum algorithms 
must be simulated on a classical machine. However, due to quantum parallelism simulating 
quantum algorithms on a classical computer comes at a potentially exponential simulation 
cost 63, 55, 40 . Therefore, as inefficiencies in simulation will only multiply rapidly dur- 
ing the evolutionary process, an efficient quantum computer simulator is a necessity in this 
field 39 . Spector 61 introduces a tested quantum computer simulator termed QGAME, 
which is based upon the quantum circuit model detailed in Section ^. 6, II Sabry outlines 
an interesting quantum simulation model in the functional language Haskell, and there are 
also many quantum computer simulators that vary in quality available on the Internet. An 
example code fragment of the quantum simulator used by Spector et al. |63| explains a possi- 
ble representation of gates in the programming language LISP, which is demonstrated below: 

(HADAMARD 0) ;; apply Hadamard transformation to qubit 

(CNOT 12);; apply CNOT gate with control (qubit 1) and target (qubit 2) 

In addition, recall that Hardy and Steeb [3T] provide some C++ and Java code for sim- 
ulating different quantum effects. Nevertheless, due to the simulating inefficiency only small 
quantum algorithms can be simulated and evolved. In this sense small quantum algorithms 
means that the following must be restricted in size: number of qubits, number of gates 
(commonly referred to as the length of a quantum circuit), and the number of iterations of 
the algorithm (with varying input) needed to determine the fitness of the algorithm. Hence, 
discovering small, yet scalable, quantum algorithms is the ultimate result from evolving 
quantum algorithms. 

4.2 Program Structures 

Within GP, the individuals (computer programs) must have a defined structure so that 
genetic operators can be applied to them. In most applications of GP, the computer programs 
consist of basic building blocks, referred to as primitives or genomes, such as constants, 
operators, problem-specific functions and inputs [31] . When evolving quantum algorithms 
the main components of the quantum algorithm are quantum gates. These primitives can 
then combine to form complete computer programs in different structures; the types of 
program structures that have been applied in the field of evolving quantum algorithms are 
tree, linear, and the hybrid linear-tree. The choice of an appropriate program structure 
appears to be a key ingredient to the success of the GP evolution |39] . 

The original GP program structure, as outlined by Koza |36l 135] . was based on the 
standard tree structure. Quantum gates and their arguments are represented as parent 
nodes and their children respectively. Children nodes can be quantum gates themselvej^l. 
and leaf nodes can be constants, zero-argument functions, qubit indices or other program 
inputs. This structure can be easily translated into a functional programming language 
where the functions are the parent nodes and the function's arguments are its children. 
The program can also be executed with a post-order traversal (left-subtree, right sub-tree, 

19 The return value of a quantum gate node is a valid index of a qubit it acted upon. For example, the CNOT 
gate is usually defined to return the index of the control qubit. 
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root), which means no further memory is required as function arguments are always locally 
accessible. It is also important to note that it is trivial to convert a tree-structured quantum 
program into the quantum circuit model, where the order of gates in the circuit model will 
be defined by the post-order traversal- 

Although Koza's tree structure provides elegant adaptations to different program sizes 
and shapes ^38^ and allows for easy expression in any functional programming language, there 
is no guarantee it is the most appropriate representation for all problems [BE]. An alternative 
to the tree structure, is a linear structure that unsurprisingly consists of a linear sequence of 
instructions. Spector et al. |66j outline reasons that suggest tree-based program structures 
offer no advantages over linear program structures, but have additional complexity. They 
state that quantum algorithms are linear, so a linear program structure appears to be more 
appropriate than a tree structure. Operators or functions in the linear based model then 
get their arguments from external memory. The linear-tree structure is a combination of 
these two structures, where there is a linear structure and branching function at each node 
of the tree structure. During execution of the tree structure, the linear structure of each 
node visited is executed. The tree structure execution begins by visiting the root node and 
executes only one path to a leaf node: the next node visited is determined by the branching 
function at the current node. In quantum computing, partial measurements of the entire 
quantum state are appropriate branching functions [39] ; therefore, the next node depends on 
the result of the partial measurement. An example branching function is if the measurement 
of qubit n yields |0) or 1 1), then the left or right subtree will be visited next respectively. 

It should also be noted that there is a graph, and consequently a linear-graph, GP program 
structure that has been applied to areas other than quantum algorithms. There are many 
types of graph structures: the well-known PADO graph structure is presented by Teller 
and Veloso |71j . It seems likely that the reason this structure has not been used to evolve 
quantum algorithms is that the complexity cost is greater than that of other structures, but 
this approach does not seem to represent quantum effects better than other structures. 

4.3 Previous Studies 

In addition to the papers cited, Leier's PhD thesis [39] is a major source for this section and 
provides more details on the specific parameters of each study than are given in this section. 
A new paper by Giraldi, Portugal and Thess [25] also contains, among other things, reviews 
of Rubinstein's [55] and Yabuki and Iba's [81] work; however, Leier's thesis [39 appears to 
be the best reference. 

4.3.1 Williams and Gray 

The book by Williams and Gray |78] is the pioneering work that suggests the use of genetic 
programming to evolve quantum algorithms. However, the goal of their study was to search 
more efficiently for alternative quantum circuits to a known quantum algorithm than con- 
ducting an exhaustive search. The difference from subsequent studies is that, in this case the 
overall unitary gate that represents the quantum algorithm was already known. This is still 
useful as it could be used to search for (hopefully faster) alternative quantum algorithms 
that solve the same problem 55 . In this case, the fitness function would include a penalty 
for circuits that were similar to the known quantum algorithm. 

Using a population size of 100, the GP algorithm implemented in this study successfully 
found a quantum teleportation algorithm (see Section [5T77T] for details on quantum telepor- 
tation), where the send portion was as efficient, and the receive portion more efficient, than 
the known algorithm; the GP algorithm used: 

• A linear program structure 

• An approximate universal-gate set consisting of the CNOT gate and two other one- 
qubit gates 

• A three-tuple representation: {quantum gate, [parameter J^Tpl set of qubits acted 
upon} 

20 There is also a trivial mapping of a quantum circuit to a tree-structured program, based on establishing a 
sequence of gates, that is further described in Leier's thesis |39j . 

21 Only discrete parameters were allowed. 

22 [ ] indicates optionality. 
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• A fitness function based on the evolved circuits similarity to the known solution 

• A ranking selection scheme 

• Mutation and Crossover operators that act on quantum gate(s) 
4.3.2 Yabuki and Iba 

Yabuki and Iba [5T] developed a GA mode0 that was specifically tailored to evolve the 
quantum teleportation algorithm as did Williams and Gray [78]. Using a larger population 
size of 5000, the GA algorithm implemented in this study successfully found a quantum 
teleportation algorithm (see Section [2T7T] for details on quantum teleportation) that has at 
least 3 fewer gates than any non-evolved quantum teleportation algorithm; the GA algorithm 
used: 

• A linear program structure 

• The same gate set as Williams and Gray [78] (chosen as the problem was the same) 

• A unique fixed-length string representation that was specifically tailored to the produc- 
tion of two entangled particles, which is essential in quantum teleportation 

• A problem specific fitness function based on the error of the evolved algorithm's output 

• A fitness-proportional selection scheme 

• Mutation and fixed-length crossover operators 

Yabuki and Iba, like Williams and Gray, have approached evolving quantum algorithms 
from an optimisation, rather than a search, standpoint. They have focused on using a 
tailored evolutionary algorithm to find a more efficient quantum algorithm for a given ex- 
isting quantum algorithm. Furthermore, along with Williams and Gray's [78] results, it has 
been shown that there is potential in this area, particularly for the quantum teleportation 
algorithm. 



4.3.3 Spector et al. 

Spector has recently published the first book [62] about using GP to automatically evolve 
quantum algorithms, which makes reference to all of his previous work with a more detailed 
introduction into the field. Spector et al. [B3] [BU [B5] [BBJ have conducted extensive research 
on evolving quantum algorithms. Three different GP models were outlined and applied to 
various problems in [66] . The common elements of all three models are the 

• Standardised fitness function (detailed in [63]) that takes into account three compo- 
nents: 

1. Hits: total number of fitness cases minus the number of fitness cases for which the 
quantum circuit yields the correct answer with less than 0.48 probability 

\ max (0. errors —0.48) 

2. Correctness: ' =1 . This formula deliberately ignores any 

maxynits , 1) " 

circuits with errors less than 0.48, so that the focus is on producing correct quan- 
tum algorithms, not on improving the probability of correctness of already correct 
quantum algorithm^] 

3. Efficiency: number of quantum gates / 100,000, which makes sure that this com- 
ponent will always be less than 1. 

Spector et al. recommend that these components be combined lexicographically, whereby 
quantum circuits are only compared on Correctness if their Hits are the same, and sim- 
ilarly only compared on Efficiency if the Hits and Correctness are the same. 

• Tournament selection method with a tournament size of five individuals 

The first GP model defined in [BBJ uses a standard tree program structure, and is in fact 
a summary of the results presented in Spector et al.'s initial paper [63] . This GP model 
was applied to two oracle problems, namely Deutsch's problem and the scaling majority-on 
problem (defined in Section |2.7.3|I . Using a population size of 10,000, this tree structure 
GP algorithm was used to evolve a better-than-classical quantum algorithm for Deutsch's 
problem, but did not evolve a quantum algorithm to solve the majority-on problem better 

It is technically not a GP model due to the fixed length representation of individuals. 

Here Spector et al. are using the term correct to mean correct greater than or equal to 48% of the time. 
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than any classical algorithm. Along with the common elements shown above, this tree 
structure GP algorithm also used 

• A set of gate building functions as well as iteration structures and arithmetic operators 

• Functions in prefix notation, analogous to functional languages such as LISP in which 
the quantum algorithms were simulated 

• Mutation operators, and two types of crossover operators: a variable length operator 
design for the tree structure and a reproduction operator to produce a child exactly 
the same as a parent. 

The other two GP models use two types of linear program structures. As stated in 
Section 14.21 a linear structure stores arguments in external memory. The first type of linear 
structure, called the stack-based linear GP, uses a global stack for temporary data storage. 
This GP model successfully evolved a quantum algorithm to solve the four-item database 
search problem faster than any classical algorithm; moreover, the evolved quantum algorithm 
is as efficient as, and almost identical to, Grover's search algorithm. Besides the program 
structure, the differences between this GP model and the tree structure model discussed 
above are that the 

• Various crossover operators are designed to operate on linear, rather than tree, struc- 
tures; they are also fixed-length operators 

• Iteration structures are stack-based, rather than tree-based 

• Gate building functions do not return values on to the stack 

• Arithmetic functions take their arguments from, and return their result to, the stack 

• A no-op operator is part of the function set 

The primary role of including iteration structures is to produce scalable quantum algo- 
rithms, however, some non-scalable quantum algorithms are of interest [66]. The second 
type of linear GP model is tailored to evolving these non-scalable quantum algorithms as 
there is no iteration structure. Furthermore, Spector et al. [66] suggest that there is no 
major reason for quantum gates to share parameter values, which consequently means there 
is no reason for data storage. Thus, in contrast to the stack-based linear structure, the sec- 
ond type of linear structure implemented, called the stack-less linear GP, only contains the 
gate building functions and has no external memory for temporary data storage. This GP 
structure was used to evolve a quantum algorithm to solve the And-Or Query problem faster 
than any previously known classical or quantum algorithm. The And-Or Query problem is 
to determine whether the boolean function (/(0) V /(l)) A (/(2) V /(3)) is true or false for 
a two-bit boolean function /. 

Spector et al. |64l 165] then presented a modification of their stack-less linear structure 
GP. The changes to the GP model include 

• Using a steady-state GP. All previous GP models applied to quantum algorithms have 
been generational GP models; the difference is simply that steady-state GP models do 
not have clearly defined generations. The remainder of the evolutionary process can 
be considered to be very similar 

• Supporting the use of variable-length representations and variable-length crossover 
operators 

• A four component fitness function, in which the components are again combined lexi- 
cographically; the components are the 

1. Misses component, which is equivalent to the Hits component used in [66] 

2. Expected- Queries component that considers the number of oracle calls (defined in 
65 ), which is a tailored modification for evolving quantum algorithms to solve 
oracle problems 

3. Max-Error component, which is the maximum probability of getting a wrong 
answer in any fitness case, and is similar to the Correctness component used in 

m 

4. Num-Gates component, which is equivalent to the Efficiency component used in 

m 

• Including a one-qubit measurement gate so that partial measurements can be made, 
which are key to several known quantum algorithms such as Shor's factoring algorithm 
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• The added ability to distribute the evolutionary process across multiple workstations 
to decrease execution time 

This modified GP model was again applied to the And-Or Query problem, with better results 
that the initial stack-less linear structure model. From these improving results, Spector 
et al. deduced that they were successfully improving their GP model to evolve quantum 
algorithms, and that the stack-less linear structure is probably the best structure that has 
been developed for evolving quantum algorithms. Furthermore, the successful research by 
Spector et al. summarised above has shown that there definitely is a degree of potential in 
evolving quantum algorithms. 

4.3.4 Rubinstein 

Using a population size of 5000, Rubinstein |55] used his generational GP algorithm to 
successfully discover the most efficient known quantum algorithms to produce two to five 
maximally entangled qubitQ, that is, qubits of the form 4= (|0 . . . 0) + 1 1 . . . 1}). The GP 
algorithm used 

• A linear program structure (with no external storage) 

• An unspecified set of quantum gates that include CNOT, Hadamard and importantly 
an 'Observe' gate that can measure one or many qubits, which is known to be a vital 
technique in several known quantum algorithms, such as Shor's factoring algorithm 

• A modification of the three-tuple representation, where the quantum gate and its 
parameters and qubit operands are encoded into a bit string, which is the standard 
representation for an individual in GPs 

• A fitness function based solely on the error of the evolved quantum system. The fitness 
was calculated using the following formula, where i is an input case of k total cases, 
j is a basis state in a quantum system of n qubits, and o and d are the observed and 
desired amplitudes respectively. 



The error obtained from this function is then divided by the highest error of any 
individual evolved to obtain a stardardised fitness that lies in the range between 
(optimal) and 1. 

• A fitness-proportional selection scheme 

• Mutation operators (with low probability), and crossover operators that act upon all 
parts of the bit string representation, that is, gates, parameters and qubits. Gate cross- 
over is variable-length, but parameter and qubit crossover are fixed-length, as most 
quantum gates act upon a fixed number of qubits with a fixed number of parameters 

From the quantum algorithms produced for the maximum entanglement problem with 
two to five qubits, deductions were made about an arbitrary sized maximum entanglement 
production circuit. Thus, this study shows that there is potential for GP to produce small, 
yet scalable, quantum circuits that can be converted into large scale quantum circuits. 

4.3.5 Lukac et al. 

Lucak and Perkowski [^3] present a general GA for evolving quantum circuits. Their algo- 
rithm used: 

• An encoding system where quantum circuits are represented as an array of strings of 
quantum gates. Each element of the array represented a specific point in time (after the 
previous element and before the next element) in the quantum algorithm when gates 
could act upon qubits. The string of ordered gates at each element corresponded to the 
order of qubits, such that the the first gate acted upon the first qubit and so on. This 
system had no extra parameters to identify which qubits a gate acted upon, thus it was 
only possible with the introduction of a no operation gate. Their encoding system also 
deliberately allowed for parallel evalution of individuals to potentially decrease fitness 
evaluation time. 

Two maximally entangled qubits is in fact the first Bell state as described in Section 12.2.51 
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• A large gate set comprising various one, two and three qubit gates, including a one- 
qubit 'wire' gate that performs 'no operation'. The reason for the diverse gate set is 
Lukac and Perkowski's focus on evolving arbitrary quantum algorithms. 

• A fitness function considering the correctness of the quantum algorithm. This fitness 
function is similar to those in previous papers 55, 78., but in this case a fitness level 
of f corresponds to maximum (not minimum) fitness. 

• Roulette wheel and stocastic universal sampling selection schemes. Stocastic universal 
sampling is less biased towards selecting 'fit' individuals, thus a more random selection 
scheme (which is detailed more in [26]). 

• Mutation and crossover operators. 

• Population sizes of 50 and f 00 individuals. The population was deliberately kept small 
to avoid long fitness evaluation times. 

Lucak and Perkowski used their algorithm to evolve various quantum circuits, averaging 
each experiment over 20 runs. First, they evolved single gate circuits (whereby the target 
gate was included in the gate set for the experiment) to test their algorithms convergence. 
All gates were successfully evolved, but for larger number of input-qubits (gate size) longer 
evolution time was required. Interestingly, mutation with a probability greater than 0.4 
was found to decrease both real-time evolution and the number of generations needed to 
successfully evolve the target gate. Their second experiment was to evolve three quantum 
circuits consisting of more than one gate: namely, the quantum teleportation algorithm 
previous evolved by Williams and Gray [78] and three and four maximally entagled qubits 
as previously evolved by Rubinstein [55]. In all three cases, the desired quantum circuit was 
evolved in similar or less time than previously published. It was also noted that while the 
number of generations required to evolve quantum circuits increases exponentially as the 
number of qubits increase, the real-time evolution increases at a much slower rate. Higher 
probability of mutation was again found to decrease the number of generations and real-time 
for successful evolution. 

While not discovering any new or further optimised quantum algorithms, Lucak and 
Perkowski have established benchmarks for the evolution of various small quantum algo- 
rithms, ranging from 1 to 4 qubits. Furthmore, for the three composite quantum algorithms 
(> I qubit) evolved, their GA has performed equally or more efficiently (in terms of time 
of evolution) than previous studies. Lukac et al. [44] have furthered their research by in- 
vestigating implementing their GA on a quantum computer. The specific encoding used in 
this GA lends itself to be computer in parallel, thus it should be implementable on quantum 
technologies, such as nuclear magnetic resonance (NMR) . 

4.3.6 Leier 

Leier, in his recent PhD thesis [39], presented two GP models for evolving quantum al- 
gorithms, one linear and one linear-tree structure. Both models are very similar to the 
successful stack- less linear structure GP model developed by Spector et al. [651 164], The 
major difference is the four component fitness function of Leier's model, in which the compo- 
nents are combined through different weightings, rather than lexicographically. The actual 
components are also slightly different, but have all been previously been mentioned by Spec- 
tor et al.; they are: misses, max-error, correctness and num-gates. Further differences in the 
linear-tree model, caused by its structural difference, are that branching functions are partial 
measurements and the crossover operators that act on the linear and tree substructures are 
included. The inherent inclusion of partial measurements and the added flexibility of the 
linear-tree structure was the motivation for creating a GP model that was not strictly linear. 

Both GP models were applied to the Deutsch-Jozsa (D-J) and f-SAT problems 26 !. The 
Deutsch-Jozsa problem is essentially a scalable version of Deutsch's problem; f-SAT also has 
a known better-than-classical quantum algorithm solution (Hogg's algorithm). The linear- 
tree GP was able to find a quantum algorithm essentially the same as the known algorithms 
for both the f-SAT and D-J problem, although some evolutionary runs did not produce a 
solution to the D-J problem. The interesting finding from this study was that intermediate 
partial measurements had no noticeable positive effect; similarly, the added flexibility of the 
linear-tree structure did not add a benefit over the strict linear structure which also evolved 

26 Note that the result of the linear-tree models applied to 1-SAT is also contained in a paper by Leier and 
Banzhaf [40]. 
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solutions to both problems. However, it is probable that these findings can not be generalised 
past these small problems that have relatively simple quantum solutions. Nevertheless, this 
study was the first to show that linear-tree structure GP models can be used to evolve 
quantum circuits. 

Leier also made some other interesting observations; the most important of these are 
listed below: 

• Using the linear structured GP model, Leier showed that there was an increase in the 
efficiency of the evolutionary process when using a pre-evolved, rather than a random, 
initial population for both the D-J and 1-SAT problem. This pre-evolution involved 
feeding evolved individuals from a smaller problem instance into the initial population. 

• According to Leier's research crossover is not as important as commonly thought. This 
finding conflicts with the traditional GP approach, where crossover is performed with 
much greater probability than mutation with the idea of multiplying and distributing 
better solutions over the population. 

• Even though Leier made a point of emphasising that his GP models were not de- 
signed to produce scalable quantum algorithms, scalability was achieved: the algo- 
rithms evolved using Leier's models, which apply to n-qubits, can be easily manually 
scaled to apply to (n+f )-qubits. However, this is in fact not surprising as the known 
quantum algorithms for both application problems are scalable. 

4.3.7 Ding et al. 

Ding et al. [T7] recently presented a new framework for evolving quantum circuits that is 
designed for both quantum algorithm discovery and optimisation. This framework uses a 
Hybrid Quantum-Inspired Evolutionary Algorithm (HQEA), which was motivated by GP 
and detailed in Yang's Masters thesis [82] . Ding et al.'s approach used: 

• A fixed length numerical encoding of quantum circuits, compared with encoding in 
symbols as done in previous works. Encoding the quantum circuits with numerical 
values was necessary to take advantage of the HQEA algorithm. 

• A gate set comprising the Hadamard, CNOT, Phase and 7r/8 gatejfj. Their approach 
is not limited to these gates, but gates in their quantum algorithms are confined to 
one-qubit and adjacent two-qubit gates. For example, the control qubit of the CNOT 
gate must be adjacent to the qubit undergoing the NOT operation. 

• A fitness function that considers both the cost and correctness of quantum algorithms 
(as also done in Reid's Masters thesis [55]). The fitness function used, where lower 
fitness is better, was: 

fitness = reward x (actualcost — satcost) + punish x (f — correctness) 

The satcost represents a satisfactory algorithms cost, whereby if it is set high or low 
then the evolution is inclined towards discovery or optimisation respectively. Algorithm 
cost was calculated with one-qubit gates, two-qubit gates, and the wire costing 1, 
2 and respectively. However, using this same framework algorithm cost could be 
more accurately computed in terms of the monetary costs using different quantum 
technologies. 

• A fitness-proportional selection scheme 

• Mutation and crossover operators. 

Ding et al. tested their approach on evolving 2 and 3 entangled qubits, as well as the 
controlled-phase gate, which confirmed that lower satcosts resulted in more optimised quan- 
tum algorithms. However, more generations of evolution were required to evolve optimised 
algorithms. Further experiments were conducted on evolving 2 entagled qubits as research 
into the appropriate values for the reward-punish factors. This revealed that a large pun- 
ish:reward ratio (> 5 for satcost = 6) is required to evolve correct quantum algorithms, and 
larger punish is required for larger satcost. In addition, Ding et al. present a faster method 
for matrix multiplication with tensor product, to increase the efficiency of evaluating indi- 
viduals. This faster method can also be adapted to other evolutionary algorithms. Overall, 
this new approach has shown promise in evolving and optimising small quantum algorithms, 
but it has not yet been used to discover a previously undiscovered quantum algorithm. 

27 Within the scope of this paper, Phase and n/8 gates can be thought of as reflection (about a basis state) and 
rotation gates respectively. 
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4.4 Further Reading 

The papers mentioned in the previous section provide further reading: Spector et al.'s work 
that was published in a book [66] is recommended as the first paper to read due to its under- 
standability. Spector's book [S2] and Leier's excellent PhD thesis |39| are recommended as 
the main points of call on the topic of evolving quantum algorithms. Both of these references 
also contain an introduction to both the field of quantum computing and evolutionary algo- 
rithms, which is similar to, but more detailed than, the one given in this paper. Spector's 
paper [61] is also recommended as an excellent short paper that extends the idea of evolving 
quantum algorithms using GP to evolving arbitrary computational processes. 

4.4.1 Related Applications of GP 

Genetic programming has also been applied to some other closely-related applications. 
Quantum program evolution based on density matrices 68 is a new area that is show- 
ing promise. Spector et al. [63] also suggested that there was potential in researching how 
GP algorithms would be executed on quantum computers, and whether a significant speed- 
up on fitness evaluation is possible using quantum parallelism. This idea has been reviewed 
by Giraldi, Portugal and Thess [25] in 2004, and Udrescu et al. |72j have presented new 
research into implementing GA algorithms on quantum computers. Another of Spector et 
al.'s suggestions is to evolve hybrid quantum-classical algorithms, whereby the quantum al- 
gorithm has classical pre and post data processors [63] . Another somewhat related paper 
[53] looked into a case of using GAs to evolve the hardware of quantum computers: a set of 
pulse sequences (which can be thought of as rotations) is evolved for a given quantum logic 
gate, implemented by Nuclear Magnetic Resonance (NMR). The significance of this is that 
the shorter and more robust the pulse sequence the more efficient the implementation of a 
quantum gate or algorithm; furthermore, the evolved set of pulse sequences was superior to 
any previously known set. 

5 Conclusion 

Genetic algorithms and programming have been successfully used to analyse and optimise 
known quantum algorithms. Previous studies have also shown that genetic programming can 
evolve new quantum algorithms, albeit only small quantum algorithms. Another possible 
outcome of evolving quantum algorithms is that a new useful idea, such as a meaningful 
sequence of gates, which will change the way future quantum algorithms are developed 
manually, may be discovered; however as yet, this has not occurred. Nevertheless, manual 
quantum algorithm generation has had more success than evolving quantum algorithms 
[39] , although no further breakthroughs have been made manually since Shor's and Grover's 
algorithms. Furthermore, the fact that quantum computers do not yet exist is a huge 
limitation on the research field of evolving quantum algorithms, as quantum simulation on a 
classical machine has an exponential order overhead. This means that only small quantum 
algorithms can be evolved, which are few in number and generally have little practical 
application. Thus, Leier and Banzhaf 39 , 40, speculate that this field will have a brighter 
future when quantum computers exist. However, there are still ideas that have not been 
tested as the field of evolving quantum algorithms is relatively new and has not yet been 
comprehensively researched |39| . 
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